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Abstract 

Measles virus (MV) causes small and large outbreaks in Iran. Molecular assays allow identifying and the sources of measles 
imported from neighboring countries. We carried out a phylogenetic analysis of measles virus circulating in Iran over the 
period 2010-2012. Specimens from suspected cases of measles were collected from different regions of Iran. Virus isolation 
was performed on urine and throat swabs. Partial nucleoprotein gene segments of MV were amplified by RT-PCR. PCR 
products of 173 samples were sequenced and analyzed. The median age of confirmed cases was 2 years. Among all 
confirmed cases, 32% had unknown vaccination status, 20% had been vaccinated, and 48% had not been vaccinated. 
Genotypes B3 and D8 (for the first time), HI and D4 were detected mainly in unvaccinated toddlers and young children. 
Genotype B3 became predominant in 2012 and was closely related to African strains. HI strains were also found in small 
and large outbreaks during 2012 but were not identical to Iranian H1-2009 strains. A majority of the Iranian D4 strains 
during 2010-2012 outbreaks were linked to the D4 strain identified in the Pakistan in 2007. We identified a single case in 
2010 belonging to D8 genotype with 99.7% identity to Indian isolates. Although the vaccination program is currently good 
enough to prevent nationwide epidemics and successfully decreased measles incidence in Iran, the fraction of protected 
individuals in the population was not high enough to prevent continuous introduction of cases from abroad. Due to 
increasing number of susceptible individuals in some areas, sustained transmission of the newly introduced viral genotype 
remains possible. 
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Introduction 

Measles is an acute and highly transmissible respiratory 
infection caused by measles virus (MV). [1] Despite effective 
vaccination that resulted in a 7 1 % drop in measles deaths between 
2000 and 20 1 1 worldwide, global mortality attributed to measles 
has been estimated to be 157700 deaths in 201 1 [2]. Human MV 
belongs to the genus Morbillivirus within the family Paramyxoviridae 
and is widespread throughout the world [1]. The C-terminal 450 
nucleotides of N-gene sequence is one of the most variable parts of 
the measles genome and commonly used for genetic character- 
ization of wild-type viruses. Based on the analysis of this genomic 
region, there are 24 recognized genotypes within the 8 clades (A— 
H) [3] . Genotyping of wild type MV allows identifying imported 
measles as well as indigenous strains [4]. Before measles 
surveillance program in 1980, measles was endemic in all 
provinces of Iran. After introduction of single-shot monovalent 
measles vaccination in 1967, a two-dose strategy (monovalent 
measles vaccine at the ages of 9 months and 15 months) was 



introduced in 1987. A mass vaccination campaign for all persons 
aged 5-25 years in 2003 (two doses of measles and rubella (MR)) 
and follow up vaccination from 2004 onwards (two doses of 
measles, mumps, and rubella (MMR) plus vaccination at the ages 
of 12 months and 18 months, measles incidence in Iran has been 
declining steadily [5]. From 2002, National Measles Laboratory of 
Iran in collaboration with Center for Disease Control (CDC-Iran) 
conducts routine molecular surveillance of MV associated with 
outbreaks and sporadic cases [5,6]. Small and large outbreaks 
have occurred periodically in Iran. Two different measles 
genotypes have been reported from 2002-May 2010 in Iran [7]. 
Genotype D4 was considered as the indigenous genotype before 
mass vaccination [8,9] and genotype HI as imported genotype 
from distant parts of Asia [7]. According to the new measles 
elimination strategy of the WHO, Iran aimed to eliminate MV by 
2015 [2]. At present, a limited number of provinces in Iran are 
prone to measles outbreaks because of inappropriate adoption of 
the measles immunization program [6] . Since the disease is highly 
contagious we aimed to obtain insight into the molecular 
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Table 1. Characteristics of the representative positive samples collected from 1 st May 2010 to 30 December 2012 in Iran*. 







Age 




Source of 




Province of 


Vaccine 


Accession 


NO Strain Designation 


(year) 


Gender 


specimen 


Genotype 


origin 


History 


number 


1 MVi/Bandarabas.lRN/5.10/1 


9M t 


F 


Urine 


D4 


South 


None 


HM440226 


4 MVi/Bandarlengeh.lRN/07.10/1 


15 


F 


Urine 


D4 


South 5 


Unknown 


HM440228 


6 MVs/Sistan.lRN/10.10 


4M 


F 


TS* 


D4 


South East 


None 


HM998697 


7 MVs/Tehran.lRN/14.10 


20 


M 


TS 


D4 


Center 11 


None 


HM998698 


8 MVs/Sistan.lRN/16.10/1 


2 


F 


TS 


D4 


South East 


Yes 


HM998699 


10 MVs/Sistan.lRN/17.10 


6M 


M 


TS 


D4 


South East 


None 


HM998701 


11 MVs/Sistan.lRN/18.10 


4 


M 


TS 


D4 


South East 


None 


HM998702 


12 MVi/Sistan.lRN/19.10/1 


4 


M 


TS 


D4 


South East 


Yes 


HM998703 


14 MVi/Sistan.lRN/20.10/1 


2 


M 


TS 


D4 


South East 


Yes 


HM998705 


16 MVi/Tehran.lRN/21.10/1 


23 


M 


Urine 


D4 


Center 


None 


HM998707 


19 MVs/lsfahan.lRN/22.10 


5 


M 


TS 


D4 


Center 


Yes 


HQ395674 


20 MVs/Sistan.lRN/23.10/1 


1 


F 


TS 


D4 


South East 


None 


HQ395675 


23 MVs/lsfahan.lRN/24.10 


3 


M 


TS 


D4 


Center 


Yes 


HQ395677 


24 MVs/lsfahan.lRN/25.10/1 


8M 


M 


TS 


D4 


Center 


None 


HQ395678 


26 MVs/Ghom.lRN/27.10/1 


16 


M 


TS 


D4 


Center 


Unknown 


HQ395680 


28 MVi/Hormozgan.lRN/28.10/1 


31 


F 


Urine 


D4 


South 


Unknown 


HQ596506 


30 MVs/Sistan.lRN/29.10 


11M 


M 


TS 


D4 


South East 


None 


HQ395681 


31 MVs/Mashhad.lRN/30.10/1 


8M 


F 


TS 


D4 


North East 


None 


HQ395682 


35 MVs/BandarJask.lRN/31.10/1 


6 


F 


TS 


D4 


South 


Yes 


HQ395684 


37 MVs/Fars.lRN/34.10 


11M 


F 


TS 


D4 


South West 


Unknown 


HQ596510 


38 MVs/Mashhad.lRN/35.10/1 


8 


F 


Urine 


D4 


North East 


Yes 


HQ596511 


41 MVs/Mashhad.lRN/39.10/1 


8M 


M 


TS 


D4 


North East 


Unknown 


HQ596514 


43 MVs/Mashhad.lRN/40.10/1 


10M 


M 


TS 


D4 


North East 


Unknown 


HQ668019 


45 MVs/Sistan.lRN/42.10 


1M 


M 


TS 


D4 


South East 1 ' 


Unknown 


HQ668021 


46 MVs/Sistan.lRN/43.10 


2 


M 


TS 


D4 


South East 1 ' 


Unknown 


HQ668022 


47 MVs/Sistan.lRN/45.10 


2 


M 


TS 


D4 


South East 


Unknown 


HQ711619 


48 MVs/Sistan.lRN/1.11 


7 


SH 


TS 


D4 


South East 


Unknown 


JF327854 


49 MVs/Fars.lRN/5.1 1 


7 


F 


TS 


D4 


South West 


None 


J F7 16447 


50 MVi/Tehran.lRN/14.11/1 


21 


M 


Urine 


D4 


Center 


Unknown 


JN048663 


52 MVs/Hormozgan.lRN/1 8.1 1 


1 


F 


Urine 


D4 


South 


None 


JN048665 


53 MVs/Sistan.lRN/52.11 


7 


F 


Urine 


D4 


South East 


Unknown 


JQ687376 


54 MVi/Yazd.lRN/2.12 


23 


M 


Urine 


D4 


Center 


Unknown 


JQ687377 


55 MVs/Tehran.lRN/9.12 


28 


M 


TS 


B3 


Center 


Unknown 


JX1 83266 


56 MVi/Sistan.lRN/12.12/1 


5 


M 


TS 


B3 


South East 


Unknown 


JX051516 


58 MVs/Sistan.lRN/14.12/1 


17 


F 


TS 


B3 


South East 


None 


JX1 83259 


60 MVi/Sistan.lRN/1 5.1 2/1 


11 


M 


Urine 


B3 


South East 


None 


JX1 83261 


65 MVi/Sistan.lRN/16.12/1 


5 


M 


TS 


B3 


South East 


Yes 


JX1 83267 


72 MVi/Sistan.lRN/17.12/1 


11M 


M 


TS 


B3 


South East 


None 


JX1 83274 


75 MVs/Tehran.lRN/18.12/1 


3 


M 


Urine 


B3 


Center 


None 


JX1 83276 


81 MVi/Sistan.lRN/19.12/1 


10 


M 


Urine 


B3 


South East 


Unknown 


JX219961 


89 MVs/Sistan.lRN/20.12/1 


5 


M 


TS 


D4 


South East 


None 


J X2 19968 


100 MVs/Tehran.lRN/21.12/1 


33 


F 


TS 


B3 


Center 


Unknown 


JX418027 


111 MVi/Sistan.lRN/22.12/1 


11M 


M 


TS 


B3 


South East 


None 


JX486009 


115 MVs/Yazd.lRN/23.12/1 


22 


M 


TS 


HI 


Center 


Unknown 


JX4 18032 


119 MVs/Ghom.lRN/24.12/1 


22 


M 


Urine 


B3 


Center 


Unknown 


JX486013 


122 MVs/Sistan.lRN/25.12/1 


14 


M 


Urine 


B3 


South East 


Yes 


JX486015 


125 MVs/Sistan.lRN/26.12/1 


10M 


M 


TS 


B3 


South East 


None 


JX631249 


127 MVs/Sistan.lRN/27.12/1 


3 


M 


TS 


B3 


South East 


None 


JX631251 


137 MVs/Sistan.lRN/28.12/1 


6M 


F 


TS 


B3 


South East 


None 


JX857307 
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Table 1. Cont. 









Age 




Source of 




Province of 


Vaccine 


Accession 


NO Strain Designation 


(year) 


Gender 


specimen 


Genotype 


origin 


History 


number 


147 MV<;/Kprman IRN/2Q 1 7/1 


1 0 


M 


75 


B3 


Cemer 


Yes 


J AOJ / J 1 Z 


i i v i v i j Lei ii - 1 ni s>\j . i z./ i 


2 


M 


jS 


B3 


•^niith Fact 

JUU LI 1 Lcl j I 




J AOJ /Jit 


148 MVs/Konarak.lRN/35.12/1 


6M 


M 


Urine 


B3 


South East 


None 


KC1 39052 


155 MVs/Chabahar.lRN/38.12/1 


2.4 


M 


TS 


B3 


South East 


Yes 


KC1 39059 


162 MVs/Chabahar.lRN/42.12 


1.2 


M 


TS 


B3 


South East 


None 


KC879273 


167 MVs/lranshahr.lRN/43.12/1 


1.2 


F 


TS 


B3 


South East 


None 


KF214765 


168 MVi/Zabol.lRN/45.12/1 


7 


M 


Urine 


D4 


East 


Yes 


KC984305 


169 MVs/Chabahar.lRN/50.12/1 


10M 


M 


TS 


B3 


South East 


None 


KC737545 


172 MVs/Konarak.lRN/51.12 


9M 


M 


TS 


B3 


South East 


None 


KC794528 


173 MVs/Zahedan.lRN/52.12 


1M 


M 


TS 


B3 


South East 


None 


KC794532 



*The nomenclature indicates the representative sample which is the first isolated strain during that week. 

1 Month. 

"'"Throat swab. 

^Pakistani who is imported cases. 
^'Afghanis who are imported cases. 
doi:1 0.1 371 /journal.pone.0094846.t001 



epidemiology of recent measles outbreaks in all provinces of Iran 
and compare virus strains to strains found in epidemics in other 
parts of the world. This information will help us to understand 
global transmission pattern of endemic measles. 

Methods 

Clinical Specimens 

The study was carried out and approved by the ethical 
committee of ministry of health of Iran for the protection of 
human subjects in Iran and written informed consent was obtained 
from all persons or their parents for participation and sample 
collection. Nose-throat swabs, urine and serum specimens of 
measles suspected cases from 1 st May 2010 to 30 th December 2012 
were collected from primary health care centers in different 
regions of Iran by Ministry of Health and Medical Education in 
Iran. All the specimens were collected 1-15 (average of 4.8) days 
after onset of rash according to WHO procedures [10] and 
transported to the National Measles Laboratory of Iran in the 
School of Public Health, Tehran University of Medical Science for 
processing. 

Antibody, RNA Detection and Viral Isolation 

All sera were screened for measles IgM specific antibody using 
the commercial ELISA kits (Dade Behring, Marburg, Germany). 
Nose-throat swabs and urine of IgM confirmed cases were also 
processed and inoculated onto Vero/hSLAM cells [1 1]. Cells were 
harvested when CPE was observed and virus isolation was 
confirmed by RT-PCR. Confirmed cases were further analyzed 
by RT-PCR and bidirectional sequencing. If virus culture was not 
successful for direct detection of MV, confirmed cases were further 
analyzed by RT-PCR and sequencing. Total RNA was extracted 
from measles isolates (infected cell lysates) or direcdy from clinical 
specimens according to the manufacturer's instructions with High 
Pure Viral Nucleic Acid Kit (Roche Diagnostics, Germany). 
Partial nucleoprotein gene segments (the C-terminus of N gene) of 
MV were amplified by RT-PCR using previously described 
primers MeV214 and MeV216 [12]. 



Sequencing and Data Analysis 

Sequencing was done on a capillary sequencer (ABI automatic 
DNA Analyzer, Life technologies) and the ABI BigDye Termina- 
tor V3.1 cycle sequencing kit (Life Technologies) according to the 
manufacturer's instructions in both directions. Wild-type measles 
isolates and genotype sequences from Iran were named as 
recommended by WHO [13]. Genotyping of the sequences was 
performed according to WHO protocol [10]. All sequence 
comparisons were based on analysis of the 450 bp sequence. 
Phylogenetic analysis using MEGA version.5 [14] was based on 
the Neighbor Joining and Kimura 2-parameter methods and 
bootstrap analyses were performed by 1,000 resampling of the 
data sets. Sequences measles strains of other countries were 
downloaded from GenBank and Measles Nucleotide Surveillance 
database (MeaNS) (http://www.who-measles.org/) and included 
in the analysis. All sequences obtained during this study were 
submitted to the GenBank and MeaNS. 

Results 

In spite of progress in measles surveillance, outbreaks were 
reported from some parts of our country over the last 3 years. A 
total of 438 cases of measles (from 1 st May 201 0-30* December 
2012) were confirmed by a serological test (IgM-ELISA). Case 
patients with measles symptoms who were IgM positive were 
considered confirmed cases. Sequence analysis was undertaken on 
173 confirmed measles cases which were both serologically positive 
and positive by RT-PCR (Table SI). 44 sequences (25.4%) were 
obtained from viral isolates using Vero/hSLAM cell culture (21 
cases isolated from throat swab and 23 from urine) and the 129 
(74.6%) were obtained directly from clinical samples. We were 
able to isolate live MV from cases during all outbreaks, which 
is useful for a more extensive characterization of the strains 
involved in the outbreaks. Demographic information for the 173 
sequences and the representative sequences are shown in the 
Table S 1 and Table 1 , respectively. The age of cases ranged from 
1 month to 60 years (6.47 ±9. 16 years). The median age was 2 
years. 37.6% (65 cases) were female and 62.4% (108 cases) were 
male. The proportion of gender was significandy different as the 
male-to-female ratio was 1.66:1. The age distribution of laboratory 
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50-i 



Unknown 
Un-vaccinated 
Vaccinated 




5-14 
Year of age 



15-29 



Figure 1. Measles, number of cases by age group and vaccination status, Iran, 2010-2012 (n = 173). 

doi:10.1371/journal.pone.0094846.g001 



confirmed cases shows 61.8% of them were less than 5 years 
(Figure 1). Age distribution of measles cases related to geographical 
location shows that the susceptible population specially in patients 
less than 5 years is increasing in South and South-East of Iran 
(Figure 2). The vaccination history showed that among all 173 
confirmed cases, 34 (20%) had been vaccinated, 83 (48%) had not 
been vaccinated, and 56 (32%) had unknown vaccination status 
(Table SI). 77% of unvaccinated cases were under 2 years old, 
who should have been vaccinated with two MMR doses. Six 
confirmed cases had traveled from Pakistan and Afghanistan 



(Table SI). Complete geographic distribution of measles cases are 
indicated in the Table SI. 69.4% of measles cases occurred in 
provinces in South-Eastern Iran that share borders with Afghani- 
stan and Pakistan. 

Phylogenetic Analysis of MV D4 Genotypes 

Sequencing results revealed that 64 (37%) out of 173 samples 
belonged to genotype D4 (Figure 3). The average sequence 
identity of 2010-2012 Iranian D4 genotypes was 94.1% and 



120 




Geographic location 



Figure 2. Age distribution of measles cases related to geographical location, Iran, 2010-2012 (n = 173). 

doi:1 0.1 371 /journal.pone.0094846.g002 
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Figure 3. Phylogenetic tree of the C-terminus of N gene 
sequences of 63selected wild-type measles isolates from Iran 
compared to the reference sequences for each genotype. The 

Iranian sequences are belonged to B3, D8, D4 and HI genotypes during 
2010-2012. Bootstrap values (1,000 replicates) >70% are indicated. 
Closely related reference sequences to the Iranian strains are shown by 
dark rhombic. Scale bar indicates base substitutions per site. Accession 
numbers of WHO reference sequences are: Genotype A (U01987), B1 
(U01998), B2 (U01994), B3.1 (AJ232203), B3.2 (L46753), C1 (AJ232203), 



D3 (U01977),D4 (U01976), D6 (L46750), D8 (AF280803), D9 (AF481485), 
D10 (AY923185), E (X84879), F(X84865), G1(U01974), G2(AF1 71 232), 
G3(AY184217), HI (AF04521 2), H2(AF04521 7). 
doi:10.1371/journal.pone.0094846.g003 

96.3% with the reference genotypes A and D4, respectively 
(Edmonston-wt.USA/54 and Montreal.CAN/89). The intra- 
sequence comparison of these strains ranged from 98.4% to 
100% (average 99.4%). One strain (MVs/Ghom.IRN/27.10/1) 
has 1.3% differences with other strains; the rest had 0.6% 
differential sequences. These data showed that the source of 
genotype D4 was not similar in all regions in 2010 because of 
multiple introductions of different D4 variants (Figure 4). All of the 
sequences were compared with the sequences previously submitted 
to GenBank and MeaNS from different parts of the world 
(Figure 4). The Iranian D4 genotypes are more closely related 
(99.3-100% identity) to wild-type MVes from Pakistan which 
were circulating in 2007 (MVs/MandiBahuddin.PAK/426.07). 
Two previous D4 strain in 2008 (MVs/Khorasan.Iran/20.08 
and MVi/Mazandaran.Iran/ 10.08) in Iran had also 100% identity 
with 2010 detected strains in Iran (MVi/Bandarabas.IRN/5.10/1, 
MVs/Sistan.IRA/10.10) and Pakistan (MVs/MandiBahuddin.- 
PAK/426.07). The results documented the continued circulation 
of genotypes D4 viruses throughout the south, east and center 
of Iran (Figure 4, 5). The MVi/Bandarabas.IRN/5. 10/2 strain 
also has 100% identity to measles isolates reported from Turkey 
(MVs/Istanbul.TUR/ 1 4. 1 1), Russia (MVs/Astrakhan.RUS/ 
52.11, MVs/Krasnodarsky-Kray.RUS/8.12), Kazakhstan (MVs/ 
South-Kazakhstan.KAZ/5.11) and Belarus (MVi/Dobrush.BLR/ 
10.12) during 201 1-2012. It raises this probability that the source of 
all were the same strain probably from Iran or Pakistan (Figure 4). 
The MVs/Ghom.IRN/27.10/1 strain has more than 99.5% 
identity with wild-type MVes from Afghanistan which were 
circulating in 2008 (MVs/Rodat.AFG/861.08). Also demo- 
graphic information showed that two Iranian measles cases 
were originally Pakistanian (MVi/Bandarlengeh.IRN/07.10/1, 
MVi/Bandarlengeh.IRN/07.10/2) and three were Afghanista- 
nian (MVs/Tehran.IRN/14.10, MVs/Sistan.IRN/42. 10, MVs/ 
Sistan.IRN/43.10) who were traveled to our country during 

2010 (Table SI). BLAST analysis revealed that the MVs/ 
IRN.Sistan/52.11 strain was 100% identical to the MVi/ 
Bandarabas.IRN/5.10/1 strain, circulating during 2010. Our 
data confirmed that susceptible persons were still present in 

2011 (Figure 4, 5). The MVi/Sistan.IRN/14.12/2 and MVs/ 
Sistan.IRN/20.12/1 are 100% identical to wild-type MVes from 
Pakistan (MVs/Islamabad.PAK/ 1 1 .08) (Figure 4, 5). 

Phylogenetic Analysis of MV D8 Genotype 

1 out of 173 samples belonged to genotype D8 (Figure 3). The 
average sequence identity of Iranian genotype D8 strain (Mvs/ 
Sistan.IRN/23.10/2) to the reference genotypes A and D8 
(Edmonston-wt.USA/54, Manchester.UNK730.94) was 94.4% 
and 97.7%, respectively. BLAST analysis of partial sequences of 
the N gene confirmed that it was 99.7% identical to the Indian 
isolates, circulating during 2010 (MVs/PUNE.IND/16.10/1). 
The patients infected with MV D8 had no history of traveling 
abroad and their transmission pathway could not be traced 
(Figure 3). 

Phylogenetic Analysis of MV B3 Genotypes 

Unfortunately the measles outbreaks intensified and continued 
to occur in 2012 in Iran. During the period from 1 st January 
through 30 th December 2012, 85 (49.13%) samples were belonged 
to genotype B3. 89.4% of these cases were belong to the south, east 
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100 



MVs/Kyzyl-Orda.KAZ/25.1 1 
MVi/Dobrush.BLR/10.12 
MVs/South-Kazakhstan.KAZ/5.1 1 
MVs/Krasnodarsky-Kray.RUS/8.12 
MVs/Astrakhan.RUS/52.1 1 
MVs/Astrakhan.RUS/1.12 
MVs/lstanbul.TUR/14.11 

♦ Mvs/lsfahan.lRN/25.10/1 

♦ Mvs/lsfahan.lRN/24.10 

♦ Mvs/lsfahan.lRW22.10 

♦ Mvi/Bandarlengeh.lRN/07.10/1 

♦ Mvi/Bandarabas.lRN/5.10/2 

♦ Mvs/Sistan.lRN/20.12/2 

♦ Mvs/Sistan.lRN/20.12/10 

♦ Mvi/Bandarabas.lRN/5.10/1 

♦ Mvs/Sistan.lRN/16.10/1 

♦ Mvs/Sistan.lRN/18.10 

♦ Mvi/Sistan.lRN/20.10/1 
♦ Mvs/Sistan.lRN/25.10/2 

♦ Mvs/IRN.Sistan/52.11 

♦ Mvi/Tehran.lRN/5.10/3 
92, ♦ Mvs/Ghom.lRN/27.10/1 

I MVs/Rodat.AFG/51.08 

♦ Mvi/Tehran.lRN/21.10/1 

♦ Mvs/Sistan.lRN/17.10 

♦ Mvi/sistan/IRN/21.10/2 

♦ Mvi/Sistan.lRN/23.10/3 

♦ Mvs/Sistan.lRN/29.10 

♦ Mvs/Mashhad.lRN/30.10/1 

♦ Mvs/Fars.lRN/34.10 

♦ Mvs/Mashhad.lRN/35.10/1 

♦ Mvs/Mashhad.lRN/39.10/1 
♦ Mvs/Sistan.lRN/23.10/1 

♦ Mvs/Mashhad.lRN/40.10/1 
♦ Mvi/Sistan.lRA/21 .10/3 

♦ Mvs /Sistan.lRN/1 0.10 

♦ Mvs/Damavand.lRN/27.12/4 
MVs/Toyama.JPN/11.12 
♦ Mvi/Sistan.lRN/17.12/2 

♦ Mvi/Tehran.lRN/14.11/1 

♦ Mvs/Hormozgan.lRN/18.11 

♦ Mvi/IRN.Yazd/2.12 
,♦ Mvi/Sistan.lRN/19.10/1 

Mvs/Sistan.lRN/19.10/2 
♦ Mvi/Zabol.lRN/45.12/1 



70 



80 



MVs/lslamabad.PAK/18.10 
MVs/Mandi Bahuddin.PAK/426.07 
MVs/Rawalpindi.PAK/6.08 

♦ Mvs/Sistan.lRN/1.11 
Mvs/Fars.lRN/5.11 

♦ Mvi/Hormozgan.lRN/27.10/2 

♦ Mvi/Hormozgan.lRN/28.10/1 

♦ Mvs/BandarJasklRN/31.10/1 
♦ Mvi/Bandarlengeh.lRN/07.10/2 

♦ Mvs/Sistan.lRN/35.10/2 

♦ Mvi/Sistan.lRN/14.12/2 

♦ Mvs/Sistan.lRN/20.12/1 
MVs/lslamabad.PAK/11.08 
MVs/Mohmand agency. PAK/21 .09 
,♦ Mvi/Sistan.lRN/16.12/2 

!♦ Mvs/Sistan.lRN/16.12/7 
♦ Mvs/Sistan.lRN/24.12/2 




— MVs/Sohar.OMN/15.12 

♦ Mvs/Sistan.lRN/42.10 

♦ Mvs/Sistan.lRN/43.10 

♦ Mvs/Sistan.lRN/45.10 
MV i/Koungou. MAY. FRA/31 .05/3 

MVi/Nairobi.KEN/23.02 
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Figure 4. Phylogenetic comparison of 51 selected Iranian D4 strains (C-terminus of N gene) during 2010-2012 with the D4 
reference sequence and other wild-type measles genotype D4 detected in other countries. Bootstrap values (1,000 replicates) >70% are 
indicated. The Iranian sequences are shown by small dark rhombic. Other wild-type measles genotype D4 detected in other countries that are not 
closely related to the Iranian D4 strains are compressed and shown by dark triangle. Scale bar indicates base substitutions per site. 
doi:10.1371/journal.pone.0094846.g004 



and south east of Iran where neighbor to the countries that have 
high measles prevalence (e.g. Pakistan and Afghanistan) (Figure 5). 
This is the first detection of genotype B3 viruses in Iran. The 
average sequence identity of Iranian genotypes B3 were 92.2%, 
98% and 96.8% with the reference genotypes A, B3.1 and B3.2, 
respectively (Edmonston-wt.USA/54, Ibadan.Nie/97/1 and New- 
York. USA/94). Although these strains belonged to the genotype 
B3, phylogenetic analysis showed that they were quite divergent 
and located in separate clusters. As demonstrated in Figure 6 
within B3 sequences at least five clusters could be distinguished. 
Iranian B3 prototypes and some other countries sequences were 
located in the cluster five (Figure 6). BLAST and phylogenetic 
analysis of partial sequences of the N gene confirmed that most of 
Iranian B3 genotype (61 strains) were 100% identical to the South 
African B3 (Cape Town) isolates, circulated during 2010 (e.g. 
MVs/Cape Town.ZAF/13.10.1) (Figure 6). According to our 
result, B3 strain (MVi/Sistan.IRN/12.12/2) was isolated form a 
Pakistani who had traveled to Iran. Some of Iranian B3 genotypes 
(9 strains) were 100% identical to the UK isolate that was reported 
at the same week in 2012 (MVs/Bradford.GBR/18.12). Also both 



are closely related (99.7% identity) to South African B3 (Cape 
Town) isolate, which circulated during 2010 (e.g. MVs/Cape 
Town.ZAF/13.10.1) (Figure 6). BLAST and phylogenetic analysis 
of the rest Iranian B3 genotype (1 1 strains) revealed that they were 
also closely related (99.7% identity) to South African B3 genotype 
(Cape Town) (Figure 6). In addition, the Iranian B3 strain i.e. 
MVs/Sistan.IRN/25.12/1 was 100% identical to the USA isolate 
that was reported at the first week in 2012 (MVs/Kansas.USA/ 
1.12/1). Both isolates are closely related (99.7% identity) to South 
African B3 (Cape Town) isolate (e.g. MVs/Cape Town.ZAF/ 
13.10.1) (Figure 6). The Intra-sequence comparison of Iranian B3 
strains ranged from 99.5% to 100% (The average was 99.85%). It 
shows that the source of many of them might be a single chain of 
transmission. 

Phylogenetic Analysis of MV H1 Genotypes 

According to the previous study, the genotype H 1 was probably 
linked to the MV importation from distant parts of Asia in 2009 
[7]. 23 (13.3%) out of 173 samples belonged to genotype HI 




Figure 5. Geographical distribution of measles virus genotypes in Iran (1 st May 2010-30 th December 2012). Numbers and type of 

genotypes are shown in boxes. 

doi:1 0.1 371 /journal.pone.0094846.g005 
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66 | ♦ Mvs/Sistan.lFW27.12/2 

14> Mvs/Bandarabas IRW29 12/2 
MVs/Ndola.ZMB/40 10/3 
MVs/Salalan.OMN/12 12 

♦ MVs/ZahedanlRN/52.12 
MVs/Alberta CANG4 10/1 
♦ Mvs/Tehran IRN/19 12/2 

♦ M\/s/Chabahar IRN/42 12 
MWSalalah OMN/17.12 

Mvs/Sistan.lRNro.^/l 

I MVs/kansas USA/1.12/1 
MV^Turkana.KEN/18.11/1 

♦ MWSistan IRW12.12/1 
M\/s/Bri1ishColumbia CAW28 12 

♦ Mvs/Sistan.lRNM 5.12/4 

MVs/Chipinge ZWE42 09 

67 i O Mvs/Hormrrozgan.lRNtte 12/5 
L- ♦ Mvs/KonaraklRN/51.12 

MVs/Mzirrba MWV23 10/2 
MVs/Matutuine MOZ/48 10 
MVs/New caslle upon Tyne GBR/26.10 

♦ Mus/Tehran.lRN'9.12 

♦ Mvs/Sistan.lRN'18.12/4 

Mi/s/Gabarone BWA/42 09 
M\/s/Rchinchua.ECU/48.11 

♦ Mvs/chabahar IRN38 12/1 
Mv/s/Tungurahua ECU/41 11/1 
MVs/Rustaq OMW06 12 
MVs/Berea LSO/51 09/1 

, ♦ Mvs/Tehran IRNM7 12/3 

G2~l ♦ Mvs/Hormozgan IRN19 12,7 
MVs/Morona ECU/10 12 
MVs/KwaZulu- Natal ZAF/24.12 
NWs/CapeTownZAF/13.1071 
MVs/Tjtume.BWA/7 10 
MVs/Minnesota USA/12 11/3 

♦ Mys/SBlan.lRN/18 12/3 

♦ Mvs/SistanlRN/20 12/3 

♦ Mvs/Sistan IRN/21 12/2 

♦ Mvs/SislanlRN/22 12/3 
MVS/Bradford GBR/18 12 



62 



81 



IWs/YoloTGO/31.10 
l*/s/Glazoue. BEN/15. 10/2 
Ml/s/Agaie.NGA/49.10/2 
MV/s/Guidan-Rouirdji NER/39 10/2 

,IWs/Nouadhibou.MRT/11.10/1 

g4 Mi/s/Bom LBR/10 10 

lWs/Grand-GedehLBR/11,10 
MVs/New York Oty.USA/18.06 
MVs/Madrid SPA/05.06/1 
Mi/s/hiarrow GBR/16 06/2 
WVs/Fbrida USA/31.06 

MVs/KhartoiimSUD/11.04 



MVs/KhartoumSULV34 02/2 
MVs/Selerrbao DR038 02/3 

MVs/Kirrtoanseke 0RO37 02/3 

MWlyon France/94 
97 1 A B3 2 {New-York USA/94) 

1 (Ibadan Me/97/1) 
I WVs/Alrreria SPA/6 03 

-MVs/MurciaSPA/18 03 

-MVs/Cela AGO/2 11 



MWNyarugenge.RWA/34.10/5 
MVs/London GBR/26 06/3 
MVs/Stutlgarl DEU/04 06 
MVi/Machakos KEN/42 05 
MVs/Freiburg DEU26 07 
MViWakisoUGA/41.06 
MViWateso.UGA/49.06 
MViWafeso.UGA/31.06/1 
MV^Texas USA/49 05 
MVs/SoesI NLD/33 07/2 
MVs/Umjiden NL034 07 

Mv/s/rJazyoona.OMN/1 5.09/1 
MVs/Mazyoona.OMN/1 6.09/1 
MVs/Mazyoona.OMN/17.09 
rA/JSalaSh.OMW19.09/1 
l*/i/Salalah.OMW19.09/2 
MUs/Sa!alab OMNT52 10 
MWSalalah.OMN/01.09 
MVs/Mazyoona OMN/1 5.09/2 




Figure 6. Phylogenetic analysis of 19 selected Iranian B3 
strains (C-terminus of N gene) during 2010-2012 with the 
B3 reference sequence and other wild-type measles genotype 
B3 detected in other countries. Bootstrap values (1,000 replicates) 
>70% are indicated. Iranian strains are indicated by dark rhombic. 
Reference strains are shown by dark triangles. Scale bar indicates base 
substitutions per site. 
doi:10.1371/journal.pone.0094846.g006 

(Figure 7). Our data showed that the average sequence identity of 
recent Iranian HI genotypes were 91.5% and 97.1% with the 
reference genotypes A and H 1 , respectively (Edmonston-wt.USA/ 
54 and Hunan.CHN/93/7). BLAST analysis indicated that recent 
Iranian HI genotypes were closely related to the Taiwan and 
China and isolates (MVs/Taoyuan.TWN/29.08 and Mvi/Shang- 
hai.PRC/ 13.07/1) (Figure 7). The Intra-sequence comparison of 



recent Iranian HI strains ranged from 99.7% to 100% (The 
average was 99.97%). It should be noted that although both 
genotypes HI detected in 2009 and 2012 were closely related 
(99.5% identity) to each other but there is no strong data to 
confirm that it was circulated from 2009 in Iran (Figure 7). 

Discussion 

According to the appropriate sample collection of measles 
suspected cases as well as serological and molecular assays, 
genotype information is always available in Iran [5,7,8]. After 
mass vaccination campaign in 2003, there were limited numbers of 
measles outbreaks in Iran from 2004 to 2008 [5,6]. However, 
measles continued to be circulated in Iran and there were several 
small and large outbreaks during 2009-2012 and most of them 
occurred in east, south and south east, specially in rural and hard 
to reach areas (CDC-Iran data). The largest outbreaks in Iran 
occurred in 2012. Genotyping showed that the most frequently 
detected genotype was genotype B3 which appeared to be the 
major strain in the center, south, east and south-east of Iran 
(Figure 5). Given the enough of data from previous studies in Iran 
it is completely obvious that the genotype B3 viruses were 
introduced by recent importation and were not endemic in Iran. 
This observation is important because the B3 genotype has never 
been previously reported in Iran. Detection of the same B3 strains 
from Iran (MVs/Sistan.IRN/18.12/3 and Mvs/Sistan.IRN/ 
25.12/1), UK (MVs/Bradford.GBR/18.12) and USA (MVs/ 
Kansas. USA/ 1.1 2/1) suggested simultaneous importations. Ge- 
notype B3 is widely distributed throughout Africa [15]. However; 
sequence identity alone is not proof for importation from African 
regions, especially when identical sequences have been detected in 
different parts of the world. 

Previous studies differentiated B3 genotype to 3 different 
clusters [15]. Based on our sequence analysis, we present 5 
clusters to differentiate measles virus within B3 (Figure 6). Our 
epidemiological data did not identify a link between Iranian cases 
with African B3 strains, indicating that how quickly MV spreads 
between countries among susceptible individuals and become 
widely disseminated in a community with high vaccination 
coverage [16]. Measles transmission can simply occur anywhere, 
and Iran is not exception of this scenario. It was previously shown 
that the genotype B3 viruses that were detected in a large outbreak 
in Africa in 2005 were soon detected throughout the United States 
and Europe [17]. There is one known case of importation from 
Pakistan with B3 genotype in week twelve 2012 at the same time of 
detecting B3 genotype in Iran, rise the probability of circulation of 
B3 in neighboring countries. Because of proper molecular 
surveillance in Iran, viruses that were indistinguishable from the 
neighbor countries were isolated during the outbreak of measles in 
Iran during 2012. 

Genotype HI was detected in 2009 [7] but there had been no 
reported cases of measles HI until 2012. Our phylogenetic analysis 
has shown that the genotype HI detected in 2012 was different 
from those detected in 2009 with 0.5% difference. It seems that 
genetic heterogeneity of the Iranian HI strains is not a result of 
increased mutation rates but is due to re-importation from South- 
East Asian countries (Figure 5,7). 

Although circulation of D4 strains have been documented 
before 2003 [8,9] and continued in the study period, the D4 
genotype was found circulated only in 7 provinces during 2010- 
2011 and 4 provinces in 2012, which indicates the limited 
circulation of D4 in the country (Figure 5). However, it cannot be 
concluded whether identical or very similar D4 viruses were 
frequently re-introduced or whether they were continuously 
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♦ Mvi/Yazd.lRN/23.12/2 
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♦ Mvs/Bandarabas.lRN/35.12/4 
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Figure 7. Phylogenetic analysis of all Iranian HI strains (C-terminus of N gene) during 2010-2012 with the HI reference sequence 
and other wild-type measles genotype HI detected in other countries. Bootstrap values (1,000 replicates) >70% are indicated. 2012- Iranian 
H1 strains are indicated by dark rhombic. 2009- Iranian H1 strains are indicated by dark circle. Reference strains are shown by dark triangles. Scale bar 
indicates base substitutions per site. 
doi:10.1371/journal.pone.0094846.g007 



circulating and cases were missed. Genetic heterogeneity of D4 
genotype in Iran might be due to multiple chains of transmission of 
this virus that is similar to circulation of the dominant D8 genotype 
in India [18]. Neighboring countries Pakistan, Afghanistan and 
Iraq have more severe problems with measles and are therefore 
likely sources of importations of D4 viruses. Sequence information 
is not regularly available in Afghanistan and Pakistan. Therefore 
the sources of genotypes are not apparent in these countries. The 
borders between Iran, Afghanistan and Pakistan are long and 
porous. Every day, some groups of people cross over these borders 
which may have brought different measles genotypes along with 
them into Iran. On the other hand, if the fraction of susceptible 
population gradually expands then the probability of a large 
measles outbreak will increase. The figures 1,2 and 7 shows that 
although vaccination program is currently good enough to prevent 
nationwide epidemics and successfully decreased measles inci- 
dence in Iran, the fraction of protected individuals in the 
population was not high enough to prevent continuous introduc- 
tion of cases from abroad. This problem was also reported from 
Central African Republic and Taiwan [19,20]. In the case of 
independent sources, our molecular characterizations alone were 
not able to differentiate between continuous circulation of 
indigenous virus and multiple introductions from the same source. 

Previously in Iran only the genotype D4 has been reported. 
However a rapid change in pattern of MV strains is observed 
when the genotypes B3 (70% of the 106 genotyped cases) and HI 
were predominant in 2012. The replacement of one prevalent 
genotype by another is a sign that vaccination coverage was high 
enough to cut off transmission of the first genotype but not high 
enough to prevent importation/establishment of another geno- 
type. This phenomenon has also been observed in other countries 
such as Ireland [21]. 

Our results have shown that 25.4% of samples were isolated 
from cell culture which is useful for a more extensive character- 
ization of the strains involved in the outbreaks. 83.2% of clinical 
samples for MV isolation were collected within less than 7 days 
after rash onset. Inability to grow virus from cell culture might be 
due to inappropriate transportation of the samples to the 
laboratory. Similar success rates were also observed in several 
studies, which were 12% [22] and 15% [23] yields. The largest 
proportion of cases occurred in toddlers (£l year) and young 
children (5-14 years). Many of the affected toddlers were younger 
than the age of the first vaccination (12 month). These 
observations highlight the importance of achieving high vaccina- 
tion coverage in both groups (Figure 1). This also raised the 
question whether the first dose of measles vaccine should be 

References 

1. Griffin DE (2001) Measles virus. Fields virology. Lippincott Williams, and 
Wilkins.1401-1441 p. 

2. Progress in global eontrol and regional elimination of measles, 2000-2011. 
(2013) Wkly Epidemiol Rec. 88(3): 29-36. 

3. Rota PA, Brown K, Mankcrtz A, Santibancz S, Shulga S, ct al. (2011) Global 
distribution of measles genotypes and measles moleeular epidemiology. J Infeet 
Dis. 204 Suppl 1: S514-523. 

4. Jayamaha J, Binns PL, Fennell M, Fcrson MJ, Newton P, ct al. (2012) 
Laboratory diagnosis, molecular characteristics, epidemiological and clinical 
features of an outbreak of measles in a low incidence population in Australia. J 
Clin Virol. 54(2): 168-173. 



brought down to 9 months again. In addition, gradually elevation 
of measles infection in vaccinated children confirms vaccination 
failure in 2-4 and 5-14 age groups in Iran (Figure 1). We believe 
that vaccination of young school aged children and toddlers can be 
necessary to prompdy control outbreaks in south and south east of 
Iran and prevent a large reservoir of susceptible persons to avoid 
new outbreaks. Age distribution of cases divided for macro-areas 
shows that the susceptible population is currendy expanding in 
South and South-East of Iran (Figure 2). In fact, a shift in the age 
distribution of cases in South and South-East of Iran could support 
our idea that the majority of cases occurred in that region is 
because of difficulties of those populations in being reached by 
vaccination as well as vaccination failure. However, in order to 
inform public health decisions how to prevent new measles 
outbreaks a computational prediction model should be developed. 
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